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Abstract 

The program Ortocartan for algebraic calculations in relativity has just been im- 
plemented in the Codemist Standard Lisp and can now be used under the Windows 
98 and Linux operating systems. The paper describes the new facilities and subpro- 
grams that have been implemented since the previous release in 1992. These are: the 
possibility to write the output as Latex input code and as Ortocartan's input code, 
the calculation of the Ellis evolution equations for the kinematic tensors of flow, the 
calculation of the curvature tensors from given (torsion-free) connection coefficients 
in a manifold of arbitrary dimension, the calculation of the lagrangian from a given 
metric by the Landau-Lifshitz method, the calculation of the Euler-Lagrange equa- 
tions from a given lagrangian (only for sets of ordinary differential equations) and 
the calculation of first integrals of sets of ordinary differential equations of second 
order (the first integrals are assumed to be polynomials of second degree in the first 
derivatives of the functions) . 



1 The motivation 

Several programs for algebraic calculations are available on the market today, and some 
of them include packages for relativity. Still, it seems that not all possibilities in applying 
such programs are known to the users and to the authors of those systems. The newest 
developments in Ortocartan include a few facilities that are not yet in general use. It is 
hoped that the present paper will demonstrate what is still possible beyond the standard 
applications. 



The implementation of the programs in Codemist Standard Lisp described in this paper was supported 
by the Polish Research Committee grant no 2 P03B 060 17. 
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2 Calculating the curvature tensors from a given met- 
ric 



This has been the main application of the program Ortocartan ever since its first release 
in 1978. It was described in several publications, most recently in Ref. 1, so it will be 
recalled only very briefly. All details can be found in Ref. 1, in the papers cited there, and 
in the latest release of the user's manual [2]. 

The main program Ortocartan takes an orthonormal tetrad of Cartan forms as its 
input data, and calculates all quantities that appear in the usual calculations in relativity: 
the determinant of this tetrad, the components of the inverse tetrad, the metric tensor 
and its inverse, the Ricci rotation coefficients, the Christoffel symbols, the Riemann, Ricci, 
Einstein and Weyl tensors, the scalar curvature. For the tensors calculated along the way it 
finds the (orthonormal) tetrad components and the coordinate components with each index 
raised or lowered as the user wishes. The user can make all kinds of substitutions, including 
those by pattern-matching (for an example of the latter see sec. 8). As an example, the 
input data for a very simple application of Ortocartan (calculating the curvature tensors 
for the spherically symmetric metric in the curvature coordinates) is shown below. 

This example, and all the other examples shown in this paper, were chosen trivial for 
illustrative purposes, so that the readers can easily see what happens and how. In actual 
research, all the programs were successfully applied also to very complicated examples. 

(ortocartan ' ( 
(SPHERICAL METRIC IN STANDARD FORM) 
(functions mu(t r) nu(t r)) 
(coordinates t r theta phi) 

(ematrix (exp nu) (exp mu) r 

(r *(sin theta))) 
(stop after ricci) 
(rmargin 60) 

)) 

(The command (rmargin 60) will fit the linelength of the output shown in the next section 
to the present text.) 

As before, the Ortocartan package includes the program Calculate that can carry out 
simple algebraic operations defined by the user. This has been described in the earlier 
publications, too, so again it is just recalled by the example below, in which it calculates 
the derivative by a; of a complicated expression that includes constants and functions. The 
input data is: 

(calculate ' ( 

(print example) 
(constants abed) 
(coordinates x) 
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(functions f (x)) 

(operation (deriv x (a ~ (b ((der x f) ~ (c ~ d)))))) 
)) 



and the output is[]: 
(print example) 

(I UNDERSTAND YOU REQUEST THE FOLLOWING EXPRESSION TO BE SIMPLIFIED) 

d 

c 

X 

b 

> deriv (x,a ) 



THE RESULT IS 

d 

c d 

f , c d 

x 

b f, d - 1 + c 

x 

> result = a b c log (a) log (b) f, f, 

1 xxx 

(I REALLY LIKED THIS ! CAN I HAVE MORE ? PLEASE ?!?) 

END OF WORK 

(RUN TIME = 50 msec) 



This is a fictitious example, not taken from any actual application, and was intended to 
demonstrate the power of our differentiating and printing procedures. 

^ach piece of output shown further is a part of a disk-output from the appropriate program that was 
inserted verbatim in the latex code of this paper. 
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3 Output in the form of Latex input and in the form 
of Ortocartan's input 

The simple call to Ortocartan shown above would result in printing all the formulae in the 
standard mathematical format. One item of the output is shown below (the tetrad i? o 
component of the Ricci tensor): 



-1 

> ricci = 2 r exp (- 2 mu) nu, + exp (- 2 mu) nu, 

r 

2 

> + exp (- 2 mu) nu, - exp (- 2 mu) mu, nu, 
r r r r 

2 

> - exp (- 2 nu) mu, - exp (- 2 nu) mu, + 
r t t t 



exp (- 2 nu) mu, nu, 
t t 



This can be either just viewed on the screen or stored on a disk and then printed. However, 
if a formula like the one above is to be inserted in a text for publication, then one may 
avoid retyping it. The same formula can be obtained in the form of Latex input code if 
one inserts the following item anywhere between the arguments of the function Ortocartan 
(but after the title (SPHERICAL METRIC...)): 



(output for latex) 

The same component of the Ricci tensor will then appear on output in the following form 



\begin{equation} 

ricci_{00} = 2r~{-l}\exp (- 2{\mu}){\nu} , _{r} + \exp (- 2{\mu}){{\nu},_{r}}~{2} 
$$ 

$$ + \exp (- 2{\mu}){\nu},_{r r} - \exp (- 2{\mu}){\mu},_{rK\nu},_{r} - \exp 
(- 2{\nu}H{\mu},_{t}}-{2} - \exp (- 2{\nu}) $$ 

$$ {\mu},_{t t} + \exp (- 2{\nu}){\mu},_{t}{\nu},_{t} 
\end{equation} 
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Note that Ortocartan has recognized the Greek letters and replaced them with the appro- 
priate Latex constructs. Then the user can add his/her favourite Latex preamble to this 
and run Latex on the result obtaining the following output: 



riccioo = 2r 1 exp(— 2/j,)u, r + exp(— 2fi)u, r 2 



+ exp(— 2ji)u. 



exp(— 2/x)//, r v, r — exp(— 2u)fi, t 2 — exp(— 2v) 



fi, tt + exp(-2u)fi, t u, t 



(1) 



The exp(f ) can be very easily replaced by using the substitutions. The exponential 
function is represented as above in order to give the user a greater freedom: the symbol 
"e" can be used for anything the user wishes and is not reserved to mean only the base of 
natural logarithms. 

If the output or some part of it is to be used as input for one of the programs of the 
Ortocartan set, then the user can avoid rewriting again. It is then enough to insert the 
following instead of the "(output for latex)" line: 

(output for input) 

The whole output will then be written in the Ortocartan input notation and any part of 
it can be inserted in the actual input just by cutting and pasting. The same component 
of the Ricci tensor that was shown above would then appear in the following form: 

ricci (0 0) = ((2 * (r ~ -1) * (exp (-2 * mu)) * (der r nu)) + ((exp (-2 * mu 
)) * ((der r nu) ~ 2)) + ((exp (-2 * mu)) * (der r r nu)) + (-1 * (exp (-2 * 
mu)) * (der r mu) * (der r nu)) + (-1 * (exp (-2 * nu)) * ((der t mu) " 2)) + 
(-1 * (exp (-2 * nu)) * (der t t mu)) + ((exp (-2 * nu)) * (der t mu) * (der t 



The automatically generated output in the input format has the tendency to use more 
parentheses than are necessary. The algebraically equivalent input written by the user 
would not contain any of the unnecessary parentheses because Ortocartan's input notation 
recognizes the usual hierarchy of priorities among the algebraic operations - see manual 
[2]. The additional parentheses do not change this hierarchy. 

The " (output for latex)" and " (output for input)" options are available for all the other 
programs described below. 



nu))) 
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4 The program Ellisevol. 



This program calculates all the quantities appearing in the evolution equations of the 
kinematical tensors of fluid flow, as defined by Ellis [3]. Since all these equations are con- 
sequences of the Ricci identity u a -^ — u a;j p = u^R^ap^, they will be fulfilled identically in 
most cases. However, they may fail to be identically fulfilled when one makes assumptions 
about separate parts of the flow, e.g. if one assumes that the shear is zero. As is well 
known, such assumptions have consequences for the other characteristics of the flow, and 
the Ellis equations will show what the consequences are. Along the way, the program cal- 
culates all those quantities that are calculated by the program Ortocartan, and in addition 
the expansion, the acceleration, the shear tensor and scalar, the rotation tensor and scalar, 
and the electric and magnetic parts of the Weyl tensor with respect to the velocity vector. 

Since the signature assumed in the calculation is (H ), the formulae may differ 

from those given in textbooks, and so some of them are quoted below for reference. The 
equations that the program verifies are the following. 

The rotation constraint equations: 



(square brackets on indices denote antisymmetrization, round brackets on indices denote 
symmetrization) . 

The shear constraint equations: 



^[a/3; 7 ] + U[ a - n U/3] + ?i[ a %] = 




where h a p = 5 a p — u a up is the projection tensor. 
The rotation evolution equations: 

ha , hp s u) y s - h^h^ii^-s] + 2a S[ a uj s p\ + -9u a p 
The Raychaudhuri equation: 



0. 




0. 



The (coordinate) electric components of the Weyl tensor: 



E a p — C ap p a u p u a = Ep, 



The shear evolution equations: 



6 S 

ho?hp <T 7 5 - ha^hp + U a Up + UJa^p + (T^CT 1 ' p + -Q<J a p 




The magnetic components of the Weyl tensor: 
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where e a7fMU is the Levi-Civita symbol. 
The "magnetic constraint" equations: 

where wp is the rotation vector field defined by: 

As an example, let us consider the application of this program to the metric of Lanczos 

[4]- 

ds 2 = (dt + Crdtp) 2 - i)d^ 2 - ^e~ r dr 2 /^ - e~ r dz 2 , 
where C = const and 

i> = (C 2 r + A-Ae- r ). 

This is a stationary cylindrically symmetric solution of Einstein's equations with a 
rotating dust source and with a nonvanishing cosmological constant A. The coordinates 
used in the metric shown above are comoving and the velocity vector field of the dust is 
one of the orthonormal tetrad vectors, hence the tetrad components of the velocity field 
are (1 0). Since this is a solution of Einstein's equations, this vector field is uniquely 
determined by the metric, and so, as expected, all the constraint and evolution equations 
will be identities. However, the acceleration (= 0), rotation, expansion (= 0), and shear 
(= 0) are all calculated, along with the electric and magnetic parts of the Weyl tensor. 

The input data is here: 



(setq !*lower nil) 
(ellisevol ' ( 

(LANCZOS METRIC) 

(coordinates t phi r z) 

(velocity 10 0) 

(constants C Lambda) 

(symbols psi = (C ~ 2 * r + Lambda - Lambda * (exp (- r))) ) 
(ematrix 1 (C * r) 

((C " 2 * r + Lambda - Lambda * (exp (- r)))~ (1 2 )) 
((12) * (exp ((-12) * r)) * (C ~ 2 * r + Lambda 
- Lambda * (exp (- r))) ~ (-1 2)) 
(exp (- (1 2) * r))) 
(substitutions (C ~ 2 * r + Lambda - Lambda * (exp (- r))) = psi ) 
(dont print messages) 
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(tensors einstein) 
)) 

(setq ! *lower t) 



The first line of the input tells the system to treat the lower-case letters and their corre- 
sponding capitals as different symbols, the last line reverses this command. The command 
"(dont print messages)" tells the system not to print the messages about unsuccessful 
attempts to carry out the substitution specified in the previous line. Without this com- 
mand, the program would write a message every time when there was no opportunity for 
this substitution in a given formula. 

The most important parts of the output are given below. Each tensor is printed with 
its unique name to facilitate substitutions. Thus the covariant rotation tensor u a/ 3 has the 
name "rotdd", while the mixed rotation tensor u a ^ has the name "rotdu", and similarly 
for shear. 



ACCELERATION = 



> rotdd = (1/2) C 
1 2 



(rotdd calculated) 
(TIME = 1060 msec) 



(rotdd completed) 
(TIME = 1070 msec) 



2 

> rotdu = - 2 C exp (r) psi 
1 



2-1 
rotdu = - (1/2) C r psi 
2 



1 -1 
rotdu = (1/2) C psi 
2 



(rotdu completed) 
(TIME = 1100 msec) 



2 

> ROTATION SQUARED = C exp (r) 
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(ROTATION SCALAR calculated) 
(TIME = 1100 msec) 



> EXPANSION SCALAR = 

(EXPANSION SCALAR calculated) 
(TIME = 1150 msec) 

(sheardd calculated) 
(TIME = 1160 msec) 

SHEAR = 

(ALL THE ROTATION CONSTRAINTS ARE FULFILLED IDENTICALLY) 
(ROTATION CONSTRAINTS calculated) 
(TIME = 1160 msec) 

(ALL THE SHEAR CONSTRAINTS ARE FULFILLED IDENTICALLY) 
(SHEAR CONSTRAINTS calculated) 
(TIME = 2500 msec) 

(ALL THE ROTATION EVOLUTION EQUATIONS ARE FULFILLED IDENTICALLY) 
(ROTATION EVOLUTION EQUATIONS calculated) 
(TIME = 2800 msec) 

> RAYCHAUDHURI EQUATION = 

(RAYCHAUDHURI EQUATION calculated) 
(TIME = 3150 msec) 

2 

> elweyl = - (1/3) C exp (r) psi 

1 1 

2 -1 

> elweyl = - (1/12) C psi 

2 2 

2 

> elweyl = (2/3) C 

3 3 

(elweyl calculated) 
(TIME = 3310 msec) 

(ALL THE SHEAR EVOLUTION EQUATIONS ARE FULFILLED IDENTICALLY) 
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(SHEAR EVOLUTION EQUATIONS calculated) 
(TIME = 3820 msec) 

> magweyl = - (1/2) C 
2 3 

(magweyl calculated) 
(TIME = 3890 msec) 

(ALL THE MAGNETIC CONSTRAINTS ARE FULFILLED IDENTICALLY) 
(magcons calculated) 
(TIME = 5170 msec) 

(I REALLY LIKED THIS! CAN I HAVE MORE ? PLEASE ?!?) 

END OF WORK 

(RUN TIME = 5170 msec) 



5 The program Curvature. 

This program calculates the curvature tensor from given connection coefficients in any 
number of dimensions. The connection coefficients are assumed symmetric (i.e. torsion- 
free), but need not be metrical. The program was written for one special application, and 
hence the assumption of zero torsion; it can be removed without any difficulty. 

6 The program Landlagr. 

This program calculates the reduced lagrangian for the Einstein equations by the Landau- 
Lifshitz method [5]. This is the Hilbert lagrangian with the divergences containing second 
derivatives of the metric already removed. In short, this (noncovariant) lagrangian is 
obtained by deleting from the scalar curvature those terms in which the Christoffel symbols 
are differentiated, and taking the remaining part with the reverse sign. 

As is well-known, the Euler-Lagrange equations derived from a variational principle 
in which the generality of the metric was limited by various assumptions (e.g. about 
symmetry) are not always equivalent to the Einstein equations. It may happen that they 
will have nothing to do with the Einstein equations; this is known, for example, to occur 
for certain Bianchi-type models [6]. Therefore, the user must make sure, when using the 
program "landlagr", that the situation is appropriate for applying the lagrangian methods. 

7 The program Eulagr. 

This program calculates the Euler-Lagrange equations starting from a lagrangian specified 
by the user. It is assumed that the resulting E-L equations will be ordinary differential 
equations (i.e. that there is only one independent variable in the lagrangian action integral). 
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As an example, the program will derive the Newtonian equations of motion for a point 
particle of mass m in the cartesian coordinates {x, y, z} from the lagrangian 

L = \n{x 2 + y 2 + z 2 ) - V(x, y, z), 

where V is a potential and x(t),y(t), z(t) are the equations of a trajectory of the particle. 
The input data is: 

(setq ! *lower nil) 
(eulagr ' ( 

(The lagrangian for the Newtonian equations of motion 
in 3 dimensions) 
(constants m) 
(parameter t) 

(functions x(t) y(t) z(t) V(x y z) ) 
(variables x y z) 

(lagrangian ((12) * m * ((der t x) ~ 2 + (der t y) ~ 2 
+ (der t z) ~ 2) - V)) 
)) 

(setq ! *lower t) 

and the results are: 

(The lagrangian for the Newtonian equations of motion in 3 dimensions) 

2 2 2 

> lagrangian = - V + (1/2) m x, + (1/2) my, + (1/2) m z, 

t t t 

(THIS IS THE VARIATIONAL DERIVATIVE BY x) 

> eulagr = m x, + V, 

t t x 

(THIS IS THE VARIATIONAL DERIVATIVE BY y) 

> eulagr = m y, + V, 

1 t t y 

(THIS IS THE VARIATIONAL DERIVATIVE BY z) 

> eulagr = m z, + V, 

2 t t z 
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(I REALLY LIKED THIS! CAN I HAVE MORE ? PLEASE ?!?) 



END OF WORK 

(RUN TIME = 100 msec) 

All the other programs described in this paper are, from the programmer's point of 
view, rather similar to each other. The various programs are simply different sets of the 
same basic algebraic operations. The program Eulagr is somewhat exceptional, since, when 
calculating variational derivatives, the functional expressions x l ,i = 1,2, ...,n have to be 
treated as independent variables at a certain stage. The program generates names for these 
variables, adds them to the list of independent variables for differentiation, replaces the 
x l by the new variables, then calculates the derivatives dL/dx l and dL/dx l , then replaces 
the new variables back by x % in the resulting expressions, and finally uses these results to 
calculate ^(dL/dx l ). All this happens automatically, and the user need not worry about 
any details of it. Ortocartan's differentiating and substituting subprograms are flexible 
enough to handle this. 

8 The program Squint. 

This program verifies whether a given expression is a first integral of a given set of ordinary 
differential equations. The program was written for a specific application, therefore it is 
rather limited in its abilities. It is assumed that the (hypothetical or actual) first integral 
is a polynomial of first or second degree in the first derivatives of the functions that should 
obey the set of equations. It is also assumed that the equations in the set are all of second 
or first order and that they have the form: 

Pitt — ^(/\ f l it , f n ,t ) 

i.e. that they are algebraically solved with respect to the highest derivatives. "Squint" is 
an abbreviation for "square integral". 

In order to make the result easy to verify, the program "squint" will be used to find 
a first integral of the equations found in the previous example. We shall at first pretend 
that we do not know what the integral I should be and will assume that it is a general 
polynomial of second degree in the first derivatives by t of the functions x(t), y(t) and z(t), 
thus / = Qij^f^f + Li^rf + E, where Qij, Li and E are unknown functions of the x l . The 
argument "(markers M)" defines M as a symbolic variable that can be used to represent 
any expression. Thanks to this facility the single equation (der t t M) = (- (der M V) 
/ m) represents the 3 equations d 2 x % /dt 2 = —dV/dx 1 for % — 1,2,3 simultaneously. The 
"maineq" directs these substitutions to the main equation, which is the explicit formula 
for the total derivative dl/dt. 

The input data is: 

(setq ! *lower nil) 
(squint ' ( 
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(a first integral of the Newtonian equations of motion) 
(constants m) 
(parameter t) 

(functions x(t) y(t) z(t) V(x y z) Qll(x y z) Q12(x y z) Q13(x y z) 

Q22(x y z) Q23(x y z) Q33(x y z) Ll(x y z) L2(x y z) L3(x y z) 

E(x y z) ) 

(variables x y z) 
(integral (Qll * (der t x) ~ 2 + 2 * Q12 * (der t x) * (der t y) 
+ 2 * Q13 * (der t x) * (der t z) + Q22 * (der t y) ~ 2 
+ 2 * Q23 * (der t y) * (der t z) + Q33 * (der t z) ~ 2 
+ LI * (der t x) + L2 * (der t y) + L3 * (der t z) + E) ) 
(markers M) 

(substitutions 
maineq 

(der t t M) = (- (der M V) / m) 
) 

(dont print maineq) 
)) 

(setq ! *lower t) 

The output includes the integral I printed in the mathematical format, the "main 
equation" dJ/dt (the printing of d//dt was suppressed by the command "(dont print 
maineq)"), the coefficients of all the expressions ^jj-^jfr^, ll~lif"> an< ^ ^he remaining 
part of dl/dt that does not contain any ^jj. This will be done correctly only if the 
second derivatives appearing in dl/dt had been replaced using the (substitutions ...), 
in this way use is made of the set of equations for which / is expected to be the first 
integral. If the set consists of n equations, then the number of coefficients printed will be 
|(n + l)(n + 2)(n + 3), i.e. 20 for n = 3. When the explicit form of / is unknown, as in 
this example, the coefficients will be partial differential equations that should determine 
the Qij, Li and E (each of the "equations" printed is a left-hand side of an equation that 
has zero on the right-hand side). 

Only a few pieces of the output are reproduced below. 

(a first integral of the Newtonian equations of motion) 
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2 

> integral = E + Qll x, +2 Q12 x, y, +2 Q13 x, z, + Q22 

t t t t t 

2 2 

> y, +2 Q23 y, z, + Q33 z, + LI x, + L2 y, + L3 z, 

t t t t t t 

> 

t 

3 

> THIS IS THE COEFFICIENT OF x, 

t 

> equation = Qll, 

1 x 

2 

> THIS IS THE COEFFICIENT OF x, y, 

t t 

> equation = 2 Q12, + Qll, 

2 x y 



[part of the output deleted] 



> THESE ARE THE TERMS THAT ARE FREE OF THE DERIVATIVES 

-1 -1 -1 

> equation =-m L1V, -m L2V, -m L3V, 

20 x y z 

(I REALLY LIKED THIS! CAN I HAVE MORE ? PLEASE ?!?) 

END OF WORK 

(RUN TIME = 1030 msec) 



Now let us substitute the well-known solution of these equations into the data and see 
what happens. 
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(setq ! *lower nil) 
(squint ' ( 

(a first integral of the Newtonian equations of motion - the final result) 
(constants m) 
(parameter t) 

(functions x(t) y(t) z(t) V(x y z)) 
(variables x y z) 
(integral ((1 2) * m * ((der t x) ~ 2 + (der t y) ~ 2 

+ (der t z) ~ 2) + V) ) 

(markers M) 

(substitutions 
maineq 

(der t t M) = (- (der M V) / m) 
) 

(dont print maineq) 
)) 

(setq ! *lower t) 
The result is: 

(a first integral of the Newtonian equations of motion - the final result) 

2 2 2 

> integral = V + (1/2) m x, + (1/2) my, + (1/2) m z, 

t t t 

(THE FIRST INTEGRAL IS ALREADY MAXIMALLY SIMPLIFIED 
AND IS EXPLICITLY CONSTANT) 

> maineq = 

(I REALLY LIKED THIS! CAN I HAVE MORE ? PLEASE ?!?) 

END OF WORK 

(RUN TIME = 130 msec) 

Similar programs for verifying first integrals can be written for any kind of dependence 
of the first integral on the first derivatives of the functions. 
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A How to obtain Ortocartan. 



The base language of this newest release of Ortocartan is Codemist Standard Lisp. The 
latter is implemented in Windows 98 and in Linux, and so should be usable for most 
computer users. 

In order to use Ortocartan one must first buy the Codemist Standard Lisp. It can be 
bought from: 

Professor John Fitch, Director 

CODEMIST Limited 

"Alta", Horsecombe Vale 

Combe Down 

BATH, Avon, BA2 5QR 

England 

phone and fax (44-1225) 837 430 
email jpffitch@maths.bath.ac.uk 

Ortocartan is free of charge. If you wish to have it, just contact A. Krasihski, I will 
either email the programs to you or send you a diskette. How to install Ortocartan when 
Lisp is already working is described in sec. 7 of the manual [2] . 
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